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Abstract 

Inelastic losses are crucial to a quantitative analysis of x-ray absorption spectra. However, current 
treatments are semi-phenomenological in nature. Here a first-principles, many-pole generalization 
of the plasmon-pole model is developed for improved calculations of inelastic losses. The method is 
based on the GW approximation for the self-energy and real space multiple scattering calculations 
of the dielectric function for a given system. The model retains the efficiency of the plasmon- 
pole model and is applicable both to periodic and aperiodic materials over a wide energy range. 
The same many-pole model is applied to extended GW calculations of the quasiparticle spectral 
function. This yields estimates of multi-electron excitation effects, e.g., the many-body amplitude 
factor Sq due to intrinsic losses. Illustrative calculations are compared with other GW calculations 
of the self-energy, the inelastic mean free path, and experimental x-ray absorption spectra. 
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INTRODUCTION 



The theory of the extended x-ray absorption fine structure (EXAFS) is now well developed 
and can be calcnlated onantitatively in m any systems. [1] However, x- r ay absorption near 
edge structure (XANES) calculations have remained only semi-quantitative at best. [2] One 
of the reasons for this disparity is a lack of accurate treatments of inelastic losses in the near 
edge region. For example, traditional calculations of EXAFS typically rely on simplified or 
semi-phenomenological models of inelastic losses in terms of a complex, energy- dependent 
exchange- correlation potential, i.e., the quasi-particle self-energy where E is the quasi- 

particle energy. In addition a many body amplitude factor Sq must be applied to the EXAFS 
signal to account for intrinsic losses, though this is frequently ignored or considered to be a 
free parameter. Two commonly used models for the self-energy in x-ray absorption spectra 
(XAS) are i) the Hedin-Lundqvist plasmon-pole model, and ii) the Dirac-Hara exchange 

n fi n n n 

approximation plus a constant complex potential. [l|, [3j, |4|, [5|, |6| Since the self-energy is 
smoothly varying at high energy and relatively small compared to the kinetic energy, these 
approximations are often adequate for EXAFS. However, variations in the self-energy tend 
to be large in the XANES region, i.e., within the first 50 eV above the Fermi energy, and 
neither of the above models describes this variation correctly. The energy scale mentioned 
above is set by the dominant excitations in the system, and is comparable to the mean 
plasma frequency u p , which is typically about 10-30 eV. Thus the EXAFS (characterized by 
weak scattering due to large loss) and XANES (characterized by large scattering and low 
loss) regions correspond to low and high energy relative to uo p . As a result, the variation 
in S(£') with energy leads to significant errors both in amplitude and peak positions in the 
XANES. 

In an effort to improve on these simplified models we present here a many-pole GW ap- 
proximation for the self-energy, 0, based on real space multiple-scattering calculations 
of the inverse dielectric function for a given system. Our goal is to develop an approach 
which can be applied routinely both for EXAFS and XANES. Analogous many-pole models 
have been used previously in calculations of the self-energy, and of the inelastic mean 



free path (IMFP), 111, [l2| with experimental optical data as input. A few first principles 



approaches that make use of pole approximations have also been developed, 
example Ref. 



15 



14 



15| For 



makes use of a band-Lanczos algorithm to calculate a many-pole approxi- 
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mation to ;he inverse dielectric matrix. For reviews of other approaches to GW calculations 



sec Refs.^ and 
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Our many-pole model yields semi- quantitative self-energies over a wide range of photo- 
electron energies from the near-edge to about 10 3 eV, which is adequate to cover both the 
XANES and EXAFS regions. 17j The approach has a number of advantages for practical 
calculations. First the method is computationally efficient in that only a few cpu hours are 
required to calculate the dielectric function, self-energy, and spectral function for a given 
system. This is significant since XANES calculations typically take several cpu-hours, while 
full GW self-energy calculations over the complete energy range of XAS experiments are cur- 
rently impractical. Finally, the approach is applicable to a wide class of materials including 
metals, insulators, and molecular systems. 

The strategy of our treatment of inelastic losses is as follows. We begin with a first 
principles calculation of the energy loss spectrum L(u) = — Im[e(qr = 0,c<j) _1 ] in the long 
wavelength limit q = 0. Q, Next this loss function is incorporated into a many-pole 
model for the self-energy which is an extension of the single plasmon-pole model of Hedin 



and Lundqvist. [71. l20l. l2lj This self-energy yields system dependent extrinsic losses due to the 



lifetime of the quasi-particle over a broad energy range. Next, to account for intrinsic losses, 
i.e., losses due to excitations of the system in response to the sudden creation of the core hole, 



we apply our many-pole model to a calculation of the quasiparticle spectral 



unction using 



an extension of the GW approximation based on the quasi-boson model. 22j, |23|, [24j This 
yields corrections to the quasi-particle approximation for XAS in terms of a convolution of 
the quasi-particle absorption spectrum with the spectral function. Moreover, the approach 
naturally includes interference terms between extrinsic and intrinsic losses and describes the 
crossover from the adiabatic- to sudden- approximation limits. 

The remainder of this paper is organized as follows. We begin with a brief description 
of the single plasmon-pole GW model for the self-energy, together with our extension to 
many poles. We then describe our approach for calculating the dielectric function at zero 
momentum transfer, as well as the extrapolation to finite momentum transfer. Next we 
compare our results for the self-energy and the IMFP with other calculations, as well as 
with experimental results for the IMFP. Subsequently we present our calculation of the 
quasiparticle spectral function and it's relation to the self-energy. We then compare our 
calculations of XAFS with experimental as well as theoretical results. Finally, we summarize 
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our results and discuss possible improvements. 



MANY-POLE SELF-ENERGY 



The many-pole model for the self-energy deve 
pole (PP) model of Hedin and Lundqvist,[7|, [20 



oped here is an extension of the plasmon- 



2l| and contains many of the same ingre- 



dients. Thus we begin with a brief description of the PP model, and subsequently describe 
the extension to a more general loss function used in this work. A more detailed description 
of the plasmon pole model is given in Appendix . Throughout this paper all quantities are 
given in Hartree atomic units (e = h = m e = 1) unless otherwise noted. We begin with the 
GW approximation for the self-energy |3] of a homogeneous electron gas in the momentum 
representation, 

d 3 q du V(q) 



E(k,E) 



E - uj — E k ^ q + i(\k — q\ — k F )5' 

where kp is the Fermi momentum. In frequency space the imaginary part of the inverse 
dielectric function (i.e., the loss function of the electron gas) is modeled as a single pole at 
u(q) = [to 2 + aq 2 + bq 4 } 1 ^ 2 , where the coefficients of the dispersion a = k F /3 and 6=1/4 are 
chosen following the prescriptions of Hedin and Lundqvist.[?j, Q This gives a single-pole 
model of the inverse dielectric function where 



- Im [e(q, uj) x ] = itu 2 p 5[lu 2 - uj(q) 2 }, 



and 



Re [e^o/p 1 ] = 1 + 



UJr 



or 



(2) 



(3) 



Inserting these results into Eq. (CQ) then yields two terms: the first term can be integrated 
analytically and gives a static Hartree- Fock exchange potential T^hf which, in the local 
density approximation (LDA) is termed Dirac-Hara exchange, 



J HF 



(k) 



kp 

7T 



1 + 



h 2 - h 2 

2kk F 



In 



kp + k 



kp — k 



(4) 



The second term, denoted by Ed(fe, E; u p ), is the dynamically screened exchange-correlation 
contribution. This contribution arises from the creation of virtual bosons. The integrals 
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over frequency and solid angle can be performed analytically (21|, leaving an expression for 
E(j(fe, E\ uj p ) in terms of a single integral over momentum transfer q. This formulation has 
been used extensively to calculate the mean self-energy E(E') within the LDA over a broad 
range of energies for EXAFS spectra. [l| The PP approximation works well at high energies 
and more generally for systems with sharp plasmon-peaks in the inverse dielectric function 
(e.g., Al), which can be described by nearly free electron gas models. On the other hand the 
model often loses accuracy at low energies for transition metals, insulators and molecules 
with more complex loss spectra, and in practice often gives unphysical structure to the 
self-energy near oj p . [l| 

In order to improve on the plasmon-pole approximation, we now introduce a more real- 
istic representation for the inverse dielectric function, using a sum over discrete poles. This 
representation preserves the analytical character of e(u;) -1 , and corresponds to a distribution 
of bosonic excitations describing the dielectric response of a material, including both inter- 
band and intraband excitations. The inclusion of this excitation spectrum in the self-energy 
naturally broadens the single PP model in a way which is characteristic of a given system. 
Moreover, the representation can be systematically improved. 

Two steps must be accomplished in order to extend the PP self-energy to a many-pole 
self-energy (MPSE): i) The first step is to obtain a suitable approximation to the energy 
loss function L(uj) = — lm[e(q = 0,co>) -1 ]. This can be done either by calculation, as is done 
here, or from experimental optical constants, ii) The second step is to extend the q = 
result to finite momentum transfer, by representing it as a weighted sum of poles, each of 
the form Eq. ([2]), which together conserve the overall strength. In addition we approximate 
the single particle Green's function G(E) as that for a free particle. This is the first term 
in the multiple-scattering expansion and ignores fine structure; hence the calculated self- 
energy represents a uniform average. With these conditions the net self-energy is simply 
the Hartree-Fock exchange contribution plus a dynamically screened exchange-correlation 
contribution which is given by a weighted sum of single pole terms, 

E d (k,E) = J2g i Z d (k,E;cu i ) (5) 

i 

with weights gi and plasma frequencies u>i. As mentioned above, E) is given by a single 
integral over momentum transfer \q\, making the calculation quite efficient. Fig. illustrates 
the self energy from our many pole model for Cu. Note that only 20 poles were needed to 
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FIG. 1: Energy loss function L(oS) = — Im[e(a;) ] in the long wavelength limit for Cu modeled 
either as a single pole (solid vertical line), or as a sum of weighted poles (vertical dashes), compared 
to the loss function as calculated by the FEFF8 code (dot dashed). 

converge this calculation, despite the relatively broad loss function of Cu. 
Inverse dielectric function 



In our approach the inverse dielectric function is calculated using the real-space Green's 
function (RSGF) method as follows: First, a modification to the RSGF code FEFF8 



25, 



26] 



is used to calculate the total absorption cross section a{ui) for a given material over a broad 
spectrum, by summing the contributions from all occupied initial states. 0, Q The results 
presented in this paper make use of atomic initial states. However current developments 
allow for the description of a continuous band of initial states within the FEFF real space MS 



framework. 



which may further improve the results. The imaginary part of the dielectric 



function e 2 is directly related to the total absorption cross section per atom cr{uS) as calculated 
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FIG. 2: Quasiparticle self-energy for Cu calculated using a loss function with fine structure (dashes) 
or without fine structure (solid). 

by the FEFF code. €2 = (n/au)a(u) where n is the atomic number density and a is 
the fine structure constant. The real part is then obtained via a Kramers-Kronig 

transform, and finally L(uj) = — Im^u;)^ 1 ] is formed by inverting e(uj). This could be a 
computationally demanding operation. However, because the self-energy involves an integral 
over e(o;) _1 , the fine structure can be neglected in all but the lowest energy part of the 
dielectric function (i.e., the first 20 eV), as shown in Fig.EJ This approximation considerably 
reduces the computational effort. It should be noted that this prescription for the calculation 
of the loss function also neglects local field effects due to the off diagonal components of the 



dielectric matrix. Nevertheless, the method 
with experiment for a variety of materials. 



las been shown to give reasonable agreement 



18 



19| Moreover, neither the self-energy nor the 



absorption spectrum (i.e. EXAFS and XANES) are highly sensitive to details of the loss 
function provided the overall weight is conserved since these quantities are given by integrals 
over the loss function. 
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Extension to finite momentum transfer 



In order to extend the inverse dielectric function to finite momentum transfer q, we 
represent the imaginary part of the loss function as a sum of closely-spaced delta functions 

L(q, w) = -Im [e{q, ^y 1 ] =7l Yl SiUJ i S ~ ^(^) 2 ) • ( 6 ) 

i 

Typically of order 10 1 — 10 2 poles are sufficient. Matching the many-pole model e(q,c<j) _1 
evaluated at zero momentum transfer to the calculation of e(c<j)~ 1 then gives the weights 
Qi and pole locations Ui respectively. Our prescription for this match is as follows: First, 
the loss function is split into N regions, each of which is represented by a single pole. For 
each region of width Aj, the pole strength and position are chosen to preserve first and 
first-inverse moments, yielding the equations defining and uji 

QiOJi = I duulm [e^)" 1 ] , (7) 

9i = --f -ImfeH" 1 ]. (8) 

For simplicity we also use the same plasmon dispersion as in the PP model. This approx- 
imation has been checked against a dispersion relation which maintains the width of the 
pole at high momentum transfers and gives similar results for materials with a broad loss 
function such as Cu, Ag, and Diamond. For materials with a sharp loss function (i.e. Al, 
Si, and Na) this approximation may not be adequate at low energies (below the plasmon 
energy) where the contribution from the particle hole continuum can dominate the loss. 

Finally, for very low energies (i.e., the first few eV where our multiple-scattering calcula- 
tions are least reliable) the calculated loss function must be corrected. For metals a Drude 
term is added and otherwise a uniform shift of the frequencies {uo{\ is carried out while 
scaling the resultant poles, so that the inverse moment matches either empirical values or 
accurate calculations of the static dielectric constant e(0), while leaving the first moment 
unchanged. 

For stability we have found it important to preserve the inverse first frequency moment, 
since this ensures cancellation of the logarithmic singularity in the derivative of Y>hf at 

n 

k = k F and E = E F .[20] This singular behavior otherwise shows up as a sharp rise in 
Re [T,(k(E), E)] within the first 10-20 eV above Ep. In metals, where cancellation is perfect, 
Re [H(k(E),E)] is fairly flat near E=Ep. In insulators, however, this singular behavior is 
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found to enhance the jump in Re [^(k^Ep), Ep)}- Thus our prescription requires a separate 
estimate either of the static dielectric function e(0) for the case of insulators or the Drude 



parameters for metals and semi-meta! 



parameterize the dielectric matrix, 



m 



s. These quantities have been used previously to 
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28 



29 



30 



31] for example Ref. 



parameters to modify a single pole model of the dielectric function, while Refs. 
generalize to a full dielectric matrix. 



271 uses simi 
and 
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ar 
29 



Our self-energy model is similar to those of Penn JJ| and Horsch et al. lOj for the valence 
contribution. One difference is that our formulation neglects the relatively small particle hole 
continuum contributions below the plasmon onset. Another is that our formulation includes 
a first order correction to the quasipartle energy as well as a renormalization constant Z 
which accounts for the quasiparticle spectral weight. Appendix gives a short discussion of 
the equivalence of the formulas for the self-energy given by Quinn (which is the starting 
point for Penn's calculations) and by Lunqvist. A further note must be made regarding the 
difference between our many-pole model and the LDA implementation in the FEFF8.2 code. 
The plasmon frequency in the current LDA model in FEFF8.2 is dependent on the electron 
density as a function of spacial coordinates. In the model discussed here, adding spacial 
dependence greatly complicates the theory and is therefore ignored. Thus our approach 
gives the spatially averaged quasiparticle correction for the whole system. We have found 
that in the XANES region, the self-energy effects on the spectrum are not sensitive to the 
density dependence. Also our calculations use the interstitial density to determine the Fermi 
momentum; the interstitial density was chosen, instead of the average density, because we 
want the model to capture the behavior of the self-energy due to interaction with the valence 
electrons. For the core electrons FEFF8 already has an option to use a non-local Dirac-Fock 



exchange which can be applied with our many-pole model. 



32] 



EXTRINSIC LOSSES 



In this section we present results which characterize the "extrinsic losses," in XAS, namely 
the self-energy and the inelastic electron mean-free-path. To confirm that our approach 
gives improved results when compared to the PP model, we have compared with other 
calculations of the self energy, including the single PP model of Hedin and Lundqvist and 
a more accurate many-pole approximation. 15|, [33J Fig. [3] shows our many-pole self-energy 
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FIG. 



3: a) Real (top) and b) imaginary (bottom) parts of the self-energy for Diamond calculated 
our many-pole model (solid); the Hedin-Lundqvist PP model (dashes); and the iterative 



using 

band-Lanczos calculation of Refs 



15l and [33j (dot dash) . 



for diamond compared with the single-pole model as well as the band-Lanczos calculation 



of Refs. 



15 and 



33. In addition, we use our results to calculate the electron inelastic mean 



free path (IMFP). 



X(E) 



(9) 



2 ME(J5)]L 

Note that this definition is not the EXAFS IMFP A E xafs,[1| 

since that quantity characterizes 
the decay of the EXAFS amplitude and includes both core hole broadening P and the 
self-energy, Aexafs = (2E)W /[\Iia£(E)}\ + T/2]. Fig. 1 shows our results for the IMFP 



for Mo, and 
algorithm 
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or comparison the single-pole model, an optical model which uses the Penn 



121 ] and experiment. 3J] Other applications of our many-pole model as well as 
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IMFP results for a number of materials have recently been presented by Sorini et al. 17] As 
can be seen in Fig. , our self-energy gives improved results for the IMFP over a broad range 
of energies. 
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FIG. 4: Inelastic mean free path for Mo calculated using our many pole model (solid); the Hedin- 
Lundqvist single pole model (dashes); a many-pole model based on optical data[ll[ [34J] (dot 
dash); and experimental data from Tanuma et al. [34| (circles) 



INTRINSIC LOSSES 

In this section we describe the treatment of intrinsic losses in a system in terms of an 
effective quasi-particle "spectral function." {22I The many-pole GW self-energy developed 
above is adequate to describe the extrinsic losses of the photoelectron in the independent 
particle (i.e., quasi-particle) approximation for the XAS. However this approximation ne- 
glects intrinsic losses due to the excitations in the absorbing medium that arise from the 
sudden creation of the core-hole. As a consequence of these excitations, the energy of the 
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absorbing photoelectron is lowered, resulting in a shift in the absorption signal. Moreover, 
one must also take into account interference between the intrinsic and extrinsic losses. Both 
intrinsic losses and interference terms can be accounted for in terms of an energy- dependent 
spectral function. 

Here we implement a many-pole model for the spectral function derived from a direct 



extension to the GW approximation and based on a quasi-boson model. 22j Within the 



approximations detailed in Ref. 



22 



the full many-body spectrum is given by a convolution 



of the single quasi-particle spectrum with an energy dependent spectral function A e s(uj , uj') , 



i.e., 



n{uS) = J dJ A cS (u,uj')fi qp (uj - lo'). (10) 

Here A e g(uj, uj') characterizes the probability density that a photon excites a photoelectron 
of energy uj — uj', as well as additional excitations (e.g., plasmons, electron-hole pairs, etc.) 
with energy uj'. 

Similarly the intrinsic many-body corrections to the EXAFS x can De represented by 
a convolution of the single quasi-particle signal Xqp an d the normalized effective spectral 
function A eS (uj, uj') (22 ] 

X(u) = / dJ A ea (uj,uj')x q p(uj - uj'). (11) 



The convolution in Eq. (ITTT) leads to a path dependent amplitude reduction in the EXAFS 
signal Sly Since the EXAFS x can De expressed as a sum of rapidly varying sinusoidal 
contributions from each photoelectron scattering path[35| with smooth amplitudes Xj{ u ) (X 
exp(2iRjk(uj)) , the amplitude reduction for each path is given by a phasor-summation 

S 2 Auj)^ [duj'A cS (uj,uj')e 2iR ^-^- k ^\ (12) 



where Rj is one half the total scattering path length of the photoelectron. As S% j(uj) is only 
weakly energy dependent, this amplitude factor can usually be approximated by a constant 
over a broad range of energies, consistent with experimental observation. In contrast the 
behavior of S 2 j(uj) for the single-pole model exhibits much more variation. 

The spectral function can be considered to be made up of a quasiparticle peak and satel- 
lites. Since broadening can be added separately, the quasiparticle part can be represented 
as a delta function of net magnitude Z e ^ at zero excitation energy, while the satellites 
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FIG. 5: EXAFS Sq{H) (upper) and the net phase <& sa t of A sat (k) (lower) for our many-pole model 
(solid) and the Hedin-Lundqvist single pole model (dashes). Note that the many-body amplitude 
and phase factors are approximately constant over a broad range of energies in the EXAFS (~ 200 
eV and above) with Sq ranging from ~ 0.91 — 0.94 and Q S at in the range —0.21 to —0.18 rad. 



represent contributions from inelastic excitations in the medium. Within the quasi-boson 
approximation 22j one has 



A cS (uj, uj') = N(v) [[1 + 2a(u)] 5(u/) + A sat (uj, u')] 



(13) 



where N(u) is a normalization constant which preserves the overall spectral weight at each 



In our approach the satellite contribution is further broken down into three terms 



corresponding to the origin of the inelastic excitation; an extrinsic part A s e ^ t coming from 
excitations created by the photoelectron, an intrinsic part Af^ arising from the excitations 
created by the sudden appearance of the core hole, and a term Aj^ er from the interference 
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FIG. 6: Satellite spectral function for a range of photoelectron momenta. From top to bottom, 
k = 16, 8, and 4 A -1 . 



between them 



(14) 



The effect of the interference tends to reduce the satellite part of the spectral function, and 
the spectral weight is shifted back to the quasiparticle peak. This variation accounts for 
the a{io) factor appearing in the weight of the quasiparticle peak. The detailed derivation 
of the components of the spectral function arising from a PP dielectric function have been 
presented elsewhere. {22] Here it is sufficient to present results characterized by the extension 
to a many-pole loss-function as in Eq. ([6]). Thus the intrinsic and interference contributions 
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are given by 



dq 



(15) 



a a;) 



2ttA; 
x In 



g^(g) 2 



x In 



Jo q^M) 2 

Ui(q) + q 2 /2 + k Q q 



(16) 



(17) 



_u>i(q) + q 2 /2 - k q_ 

where k = [2(u — a/)] 1 / 2 is the photoelectron wavenumber and k = [^(u;)] 1 / 2 is the on-shell 
photoelectron wavenumber. The extrinsic contribution to the spectral function can be found 
from the photoelectron self-energy E(/c, cu + u/) and renormalization constant Z&. 



4£M 



n\Z, 



ko I 



r fc + ImS(A;,a; + cj / ) 



Im Z, 



feo e -(^'/2^) 2 



a;' 



(18) 



where 



A fc = Re[E(ifeo,w) - Z(k,u + lu')] 
Tk = —lm[H(ko,u)—H(k,u} + u}')]. 



(19) 
(20) 



XAS CALCULATIONS 



As illustrations of our approach, we now compare our results for the XANES spectra 
of Cu and diamond as calculated with the many-pole model against those calculated with 
the single-pole model, and with experiment. Calculations of the self-energy and the spec- 
tral function were converged with respect to the number and distribution of poles used to 
represent the dielectric function. We find that typically only 10 — 20 poles are needed to 
represent a relatively broad loss function such as that for Cu. The full multiple-scattering 
FEFF8 calculation for Cu was converged with respect to the cluster size as well as the angu- 
lar momentum cutoff l max - There are only two free parameters in our calculations; a small 
imaginary shift in the potential was used to account for experimental broadening; and a real 
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energy shift was introduced to correct for the inaccuracy in Fermi energies calculated by the 



FEFF8 code. 



36] 



Comparison with experiment in Fig. [7J and M shows a clear improvement both in the phases 
and amplitudes of the XAFS signal and the near-edge structure. These improvements can 
be linked respectively, to the real and imaginary parts of the self-energy. The real part 
induces phase shifts in the signal while the imaginary part is directly related to the inelastic 
mean free path and hence the amplitudes. Fig. [7J shows our Cu K edge XANES calculations 
with both single and many-pole self-energies compared to experiment. A large (500 atom) 
cluster was used to calculate the spectra up to ~ 35 eV above the Fermi level, above which a 
smaller (177 atom) cluster was used with higher angular momentum components. Thus we 
used l max = 3 for low energies and l max = 4 above ~ 35 eV. This was done in order to ensure 
that errors due to finite cluster size and angular momentum cutoff were small compared to 
effects of the self-energy on the XANES spectrum. The result shows improvement in the 
amplitudes and phases of the peaks, especially in the region from « 10 — 50 eV (top). The 
amplitude of the "whiteline" peak (a) is substantially reduced by the corrected self-energy, 
while the second peak (b) acquires a phase shift. The dip seen at m 32 eV (c) also attains 
a significant phase shift and an increase in amplitude. The considerable improvements 
seen in this low energy XANES region can be attributed to the fact that the plasmon pole 
self-energy has singular behavior near the plasma frequency. This behavior is absent in 
the many-pole self-energy which is naturally broadened by the width of the loss function. 
Furthermore, there is improvement even in the EXAFS region 45 — 80 eV (bottom). Here 
the single plasmon pole model gives a smooth, almost featureless curve, whereas the many- 
pole model as well as the experiment show noticeable features at 62 eV (d) and 71 eV 



(e). Fig. [8] presents similar calculations for t 
from non-resonant inelastic x-ray scattering 



re diamond K edge XANES compared to data 
38 |. For diamond we could not fully converge 



the multiple scattering calculations with respect to cluster size at all energies because of 
memory requirements of the code. Thus we present our results for a 500 atom cluster with 
lmax = 2. Here the results are more difficult to interpret because of errors due to finite cluster 
size and our approximate treatment of core hole effects. However, qualitative improvement 
is seen in the amplitudes of the EXAFS from m 25 eV on. Specifically, the feature seen in 
the experiment at approximately 32 eV (a) is absent in the single plasmon pole calculation, 
but appears in the new calculation. Also, the three subsequent peaks (b, c, and d) are 
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FIG. 7: a) Top: Cu K edge XANES calculated from the many-pole self-energy and spectral 
function of this work (solid), and for comparison the conventional single-pole model (dashes), and 



experiment 



3J (+); b) Bottom: Cu K edge EXAFS shown from k = 3.75 - 4.55A -1 . 



enhanced as a result of the new many-pole calculation, giving better qualitative agreement 
with experiment. To reiterate, the single plasmon pole self-energy has a sharp turn-on of the 
imaginary part which saturates too early, giving excessive broadening in the range beyond 
the plasmon energy. Similar self-energy effects have been seen in the F K-edge spectrum of 
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LiF, where a more computationally demanding full-GW calculation was performed. 



33J 




FIG. 8: Diamond K edge XANES calculated with the many-pole self-energy and spectral function of 
this work (solid), and for comparison the conventional single-pole model (dashes), and experiment 
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■ m 



In addition to our XANES calculations, we have performed a comparison of 



experiment 



and theory of the Cu K edge EXAFS using the analysis software ATHENA. [39| To reduce 
Debye Waller effects, which are highly correlated with the effects of the self-energy and 
many-body amplitude reduction factor, we used data taken at a low temperature of 10 K. 
In order to give a fair comparison of the two theories (PP self-energy and MPSE) we have 
fixed all parameters to empirical or theoretical values. First the theory and experiment were 
aligned by matching features in the range — 300 eV. Then background subtraction and 
normalization was performed using the same spline fitting range as well as normalization 
range for experiment and theory. The EXAFS x{k) was Fourier transformed using a k- 
range of 2.632 A -1 to 15.5 A -1 with a weighting of k 1 . Debye waller factors were set using 
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FIG. 9: Comparison of Cu K edge EXAFS Fourier transform x(i?). Theoretical calculations using 



3J (+). The 



the many-pole self-energy (solid), the plasmon-pole model (dashed), and Experiment 
value of Sq was found to be ~ 0.93 over most of the EXAFS range. 



the correlated Debye model with Od = 315K. In addition, the theory was broadened by 
0.45 eV half width half max to account for experimental broadening. The estimate of the 
experimental resolution was obtained by comparing to the width of the edge step. As can be 
seen in Fig. |9]the amplitude of the first shell peak is reduced by our new treatment of inelastic 
losses, thus improving the agreement with experiment and demonstrating the adequacy of 
our calculation of Sq. Our value for Sq (~ 0.93) also agrees with a crude approximation 
(previously implemented in the FEFF8 code) which calculates the many-body overlap of the 
atomic system and gives Sq = 0.95. 
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CONCLUSIONS 



An efficient many-pole model for calculations of inelastic losses has been developed and 
successfully implemented in an extension of the multiple scattering code FEFF8. Our many- 
pole model is based on an ab initio calculation of the zero momentum transfer loss function 
by means of the RSGF approach implemented in an extension of the FEFF8 code. Extrapo- 
lation to finite momentum transfer is performed by representing e(q, uj)" 1 as a sum of poles. 
The approach yields both the quasiparticle self-energy to account for extrinsic losses and 
the many-body amplitude factor Sq to account for intrinsic losses and interference terms. 



The validity of the sel 
principles calculations. 



o 

energy model was checked by comparison with more detailed, first 



]j| We find that Sq is nearly energy independent over a broad range. 
Calculations with the many-pole model are shown to improve agreement with experimental 
results for the near edge XAS of several materials. In addition we find that our model 
is consistent with the plasmon-pole model when applied to the extended (EXAFS) region, 
which is an important step toward quantitative full spectrum calculations. A drawback of 
the present model is that it does not fully account for the contribution due to the particle- 
hole continuum at low energies. Thus the current approach may be expected to give better 
results for materials with broad loss functions, since in these cases the self-energy will be 
dominated by plasmon-like excitations even at low energies. Other improvements would be 
to represent the energy loss function as a sum of broadened poles with momentum trans- 
fer dependent broadening, and to better account for the particle-hole continuum; efforts 
along these lines are in progress. Finally we note that the development of ab initio calcu- 
lations of inelastic losses here, together with the recently developed ab initio Debye Waller 
factor calculations 40J yields improved first principles calculations of XAS from structural 



coordinates alone, without phenomenological models or the need to fit theoretical model 
self-energy, mean free path, or many-body amplitude parameters. 

This work was supported by DOE Grant DE-FG03-97ER45623 (JJR,MPP), NIH NCRR 
BTP Grant RR-01209 (JJK), NIST Grant 70 NAMB 2H003 (APS), and Academy of Finland 
Contract No 201291/205967/110571 (JAS) and was facilitated by the DOE Computational 
Materials Science Network. 
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PLASMON-POLE SELF-ENERGY 



Here we give a more complete description of the plasmon-pole (PP) model of Hedin and 



Lundqvist.[7, 
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2l| We begin with the GW approximation for the self-energy [3] 

/d 
— G(r,r',E-uj)W(r,r',Lu). 
27T 



(21) 



Here G is the single-particle Greens function, which has a spectral representation 

G(r,r,E)-^ E _ Ei + i6sga{Ei _ EF) , (22) 
and W is the dynamically screened coulomb potential, 



W(r,r',w) = J d 3 r"e(r,r",^)-V(r",r'), 
V(r,r') ' 



(23) 

, - , . (24) 
\r-r'\ y 1 

Here V is the bare Coulomb potential and e _1 is the inverse dielectric matrix. Using the 
Green's function for a homogeneous electron gas, the self-energy in the momentum rep- 
resentation is given by Eq. (JT]). In frequency space the imaginary part of the inverse di- 
electric function (i.e., the loss function of the electron gas) is modeled as a single pole at 



u(q) 



+ aq 2 + 6g 4 ] 1//2 , where the coefficients of the dispersion a = kp/3 and 6 = 1/4 



are chosen to give the Thomas- 

n n 



ermi potential at low frequency, as well as the correct high 



momentum transfer limit. [7|, |2l| This gives an inverse dielectric function whose imaginary 
part is given by 

L(q, u) = -Im [e(q, u)- 1 } = nu 2 p 5[uj 2 - u(q) 2 }. (25) 
The real part of the loss function can be obtained via a Kramers-Kronig transform 



Re [e^r 1 ] 



7T 



duj' 



2u' 



LO 2 — oo' 2 



Im [e^,^')" 1 ] 



(26) 



lu 2 — tu(q) 2 

Inserting these results into Eq. (CQ) then yields two terms: the first term can be integrated 



analytically and gives a static Hartree-Fock exchange potential 



^HF(k) = - — 
7T 



1 + 



2kk f 



■In 



kp + k 



k F — k 



(27) 
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The second term, denoted by £<i(fc, E; u p ), is the dynamically screened exchange-correlation 
contribution, which can be interpreted as the one loop diagram containing the electron 
propagator G and a boson (plasmon) propagator 

D to>»)= 2 2 mL-, - (28) 

Thus the dynamic term arises from the creation of virtual bosons which interact with 
the photoelectron via an effective coupling |g(<?)| 2 = uj 2 V(q)/2uj(q). The integral over u 
and solid angle in Eq. ([1]) can be done analytically, so that E) is given by a one- 

dimensional integral over \q\. The resulting expression for the self-energy is quite lengthy 
and is not reproduced here, but can be found in Eq. (13) of Ref. [21 



EQUIVALENCE OF SELF-ENERGY FORMULAE 

In this appendix we demonstrate that self-energy expressions of Hedin and Lundqvist and 
of Quinn and Farrell are essentially equivalent, except for slight differences the approxima- 
tions used. We start with the self-energy of an electron gas within the GW approximation as 



given by Hedin and Lunqvist [20J] in Eq. (TO . This expression can be split into two terms. The 
Hartree-Fock exchange potential T,hf, and the dynamically screened exchange-correlation 
potential which includes the dynamic response proportional to [e(qr,cij) _1 — 1]. 

Z(k,E) = X HF (k)+X d (k,E) (29) 

and 

d 3 q duj 

v 

1 



^,£)=i/(g 5 ^W[e(9,-)- 1 -l] 



E — oj — E k _ q + i(\k — q\ — k F )5 
If we rewrite e{q 1 uj)~ 1 in its spectral representation 

2c/ 



x 



2vr oj 2 - (a/ - i5) 2 E — uj — E k . q + ifljfe - q\ - k F )5' 
22 



(30) 



/Iuj 
duj' - 2 Im [eiq^y 1 ] , (31) 

Eq. (|3"Uj) becomes 

E d (k, E) = - l -£ j^L 2u'V{q)\m [e(q, u')' 1 ~ l] 



(32) 



The integral over uj can be performed by deforming the contour to the imaginary axis and 
including residues of the Greens function when necessary. The integral along the imaginary 
axis is purely real, thus the imaginary part of the self-energy is given by the imaginary part 
of the contribution from the residues of the poles in the Greens function. The result can then 
be split into two terms. One arises from the particle contribution and occurs for energies 
greater than the Fermi energy 

d 3 q 



Im [E(k,E)} 



(2tt> 



r 0(lm [AE k _ q })Q(Re[AE k _ q }) 



x — Im [e(q, AEk.q)- 1 ] 



(33) 



(34) 



where 

AE k _ q = E- E k _ q + i6(\k - q\ 

The other is associated with the hole contribution where the energy is less than the Fermi 
energy 

d 3 q 



Im [E(k,E)\ 



(-Im [AEk-q]) (-Re [AE k _ q ]) 



1 



(2tt) 3 
Im [e(q, AE^)- 1 ] 



(35) 



Quinn and Ferrell make the further approximation that E = fc 2 /2, which yields the formula 



derived in Ref. 



41 



and is used as a starting point 



is equivalent to that of Hedin and Lundqvist 
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w Penn. l_l|, |42[ Thus, Penn's formulation 



2l| with zeroth order approximations for 



the quasiparticle energy and renormalization constant. 
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